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ABSTRACT 

The constituents of the cosmic IR background (CIB) are studied at its peak wavelengths (100 and 160 pm) by exploiting 
Herschel/PACS observations of the GOODS-N, Lockman Hole, and COSMOS fields in the PACS Evolutionary Probe (PEP) 
guaranteed-time survey. The GOODS-N data reach 3<r depths of ~3.0 mjy at 100 pm and ~5.7 mjy at 160 (im. At these levels, 
source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160yum. Differential number counts ex- 
tend from a few mJy up to 100-200 mjy, and are approximated as a double power law, with the break lying between 5 and 10 mjy. 
The available ancillary information allows us to split number counts into redshift bins. At z < 0.5 we isolate a class of luminous 
sources (Li R ~ 10 11 L Q ), whose SEDs resemble late-spiral galaxies, peaking at ~130/im restframe and significantly colder than what 
is expected on the basis of pre-Herschel models. By integrating number counts over the whole covered flux range, we obtain a surface 
brightness of 6.36 ± 1.67 and 6.58 ± 1.62 [nW nT 2 sr" 1 ] at 100 and 160 fan, resolving -45% and -52% of the CIB, respectively. 
When stacking 24 pm sources, the inferred CIB lies within l.lcr and 0.5cr from direct measurements in the two bands, and fractions 
increase to 50% and 75%. Most of this resolved CIB fraction was radiated at z < 1.0, with 160 pm sources found at higher redshift 
than 100 pm ones. 
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The c osmic IR background (CIB, Puge t et al.ll996trHauser et ail 

1998) accounts for roughly half of the total extragalactic back - 

! ground light (EBL. lHauser & Dwekll200ll:lLagache et alJl2005l) . 

half of the energy radi ated by all galax ies, at all cosmic 



i.e. 



epochs, at any wavelength (Dole et alJ 1200 6). It is therefore a 
crucial constraint on modes and times of galaxy formation. 

Deep cosmological surveys carried out with t he In frared 
Space Observatory (ISO, see Genzel & Cesarsky 2000 for a 



summary) and Spitzer Space Telescope dSoifer et alj|2008l for 
a review) produced large samples of mid-IR sources a nd deep 
number counts (El baz et all 2002; Papovic h et alj|2004l) . These 
surveys led to mid-IR CIB lower limits within a factor of two 
from the upper consttaints set by TeV cosmic opacity measure- 
ments (e.g. France schini et al . 2009). 

However, at CIB peak wavelengths (100 - 200 pm), the na- 
ture of individual galaxies building up the EBL is poorly known. 
Past surv eys produced limited samples of distant far-IR ob- 
jects (e.g. lFraver et al.l|2009l) , mainly due to the small apertures 
of the available instruments and the low sensitivity in the far- 
IR. In the 160 pm Spitzer/MIPS band, -7 % of the CIB wa s 
resolved into individually detected objects dDole et alJ 12004). 
and it was only through stacking the 24 pm source s that most 
(60 - 70%) of the far-IR CIB could be recovered dDole et alJ 
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2006; Bet hermin et alJ l2010). Similarly, at longer wavelengths, 
only stacking of 24 pm sources on BLAST maps could ac- 
coun t for the m ajority of the EBL at 250, 350, and 500 pm 
( Ma rsden et aT]l2009l) . 

With the favorable diffraction limit of the large Herschel 
3.5 m mirror (Pilbratt et al., 2010), and the high sensitivity of 
its Photodetector Array Ca mera (s r Spectrometer (PACS, 70, 
100, 160 pm; Poglitsch et al.. l2010h . confusion and blending of 
sources are much less of a limitation. We are now able to resolve 
a large fraction of the CIB at its peak into individual galaxies. 

The PACS Evolutionary Probe (PEP) extragalactic survey 
samples 4 different tiers: from the wide and shallow COSMOS 
field, through medium size areas like the Lockman Hole, all the 
way down to the 160 and 100 pm confusion limit in the pencil- 
beam, very deep observations in GOODS-N and GOODS-S, and 
even beyond by exploiting gravitational lensing in low-redshift 
galaxy clusters (e.g. Abell 2218, Altieri et al. l2010l) . 

Here we exploit the science demonstration phase (SDP) 
observations of the GOODS-N field, complemented with the 
COSMOS and Lockman Hole wider and shallower layers of the 
survey, to build far-IR galaxy number counts and derive the frac- 
tion of CIB resolved by PEP. We then take advantage of the ex- 
tensive multi-wavelength coverage of the GOODS-N field, with 
the aim of identifying the CIB contributors and the epoch when 
the bulk of the Universe IR energy budget was emitted. 
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Fig. 1. Number counts at 100 {left) and 160 fim (right), normaliz ed to the Euclidean slope. Filled/open symbols belong to flux 
bins above/below the 80% com pleteness limit. Models belong to Lagache et al1 (l2004l) . lFranceschini et al] ( l2009l) . lRowan-Robinsonl 
(120091) . iLe Borgne et al.l (12009 ), Valiante et al. (2009), L acey et alJd200S ), Gruppioni & Pozzi (in pre p.). Shaded areas represent 
ISO and Spitzer data (jRodighiero & Franceschini 2004; Herau deau et al.ll2004t IBethermin et al.ll2010l) . Inset: Collection of PACS 
number counts, including this work (red, limited to rela tive errors <15% for clarity sake), PEP Abell 2218 (grey, Altieri et al. 
i and HerMES-PACS (light-blue, Aussel et al. l2.Q10l) . The solid lines in the insets mark the trends expected for a non-evolving 



population of galaxies. 



2. Number counts 

The area used to derive number counts includes observations 
in the GOODS-N (-140 arcmin 2 ), COSMOS (-2 deg 2 ) and 
Lockman Hole (LH, -450 arcmin 2 ) fields. Appendix [A] de- 
scribes data reduction and catalog extraction. 

PACS number counts were computed in each field within 0. 1 
dex flux bins, and are reported in Fig.Q~]and TableQ] normalized 
to the Euclidean slope (dN/dS oc S~ 25 ). Incompleteness and 
spurious sources have been dealt with using Monte Carlo sim- 
ulations (see Appendix lAl. Filled symbols belong to the con- 
servative 80% comple t eness limit. Thanks to the method de- 
scribed in ICharv et all d2004j) . we extend the analysis down to 
the faintest objects detected at the 3<x level. Error bars include 
Poisson statistics, flux calibration uncertainties, and photomet- 
ric errors. The latter have been propagated into number counts 
via 10 4 realizations of random Gaussian flux errors applied to 
each PACS source, using a dispersion equal to the local mea- 
sured noise. The 95 fim ISO counts dRodighiero & Franceschini! 
2004t iHeraudeau et al.l|2004l) and the 160 fim Spitzer results by 



Bethermi n et al.l (120101 including GOODS/FIDEL, COSMOS 
and SWIRE fields) are consistent with PACS data, excepted the 
ISO faintest flux bin. 

Small insets in Fig. Q] show results from PACS Guaranteed- 
Time surveys: this work, Abell 2218 lensed counts (Altieri et 
al. l2010t) . and HerMES (Aussel et al. l2010l) . The three datasets 
complement each other and show overall agreement, the 160 /mi 
HerMES counts excepted. This discrepancy can be ascribed to 
the sheet of galaxies known to bias the HerMES LH North field, 
on which their counts are based ( Owen & Morris on 2009). 

Differential, Euclidean-normalized counts display various 
well-known features: below -100 mJy at 100 fim and -200 mJy 
at 160 fim, the slope is super- Euclidean, reaching a peak in the 
counts between -5-10 and 30 mJy, while at fainter fluxes 
the slope quickly becomes sub-Euclidean. These attributes are 
generally interpreted as indicating evolving properties in the 
IR galaxy population (see Sect. 0): the trend expected for no- 
evolution is shown in Fig. Q] insets. 



The non-normalized differential counts can be represented 
by a power law of the type dN/dS oc S a , with two distinct 
slopes at fluxes fainter/brighter than S break- The results of a 
weighted least squares fit are given in Table [2] for the different 
PEP fields. The S break flux is -5.0 mJy at 100 fim and -8.5 
mJy at 160 //m. At the bright end, the number of sources in 
the small fields is low, and error bars are dominated by Poisson 
statistics. Consequently, the large uncertainties allow for nearly- 
Euclidean fits in GOODS-N and the Lockman Hole. In the case 
of COSMOS, the slope is significantly steeper than Euclidean at 
160 fim, while the 100 fim counts are still very flat. At the faint 
end, the 160 fim slope derived from our counts (GOODS-N only) 
is consistent with wh at is found by stacking 24 fim sources on 
Spitzer 160 fim maps (IBethermin et al.ll2010i a = -1.61 ±0.21). 



Field 


Flux range 


Slope 


Error 


band 


[mJy] 


a 


da 


GOODS-N 100 


2.8-5.6 


-1.32 


±0.27 


GOODS-N 100 


5.6-45 


-2.31 


±0.05 


LH100 


6.0-89 


-2.60 


±0.05 


COSMOS 100 


8.0-142 


-2.63 


±0.03 


GOODS-N 160 


5.6-9.0 


-1.67 


±0.32 


GOODS-N 160 


9.0-71 


-2.54 


±0.04 


LH160 


9.0-113 


-2.68 


±0.05 


COSMOS 160 


17.0-179 


-3.02 


±0.04 



Table 2. Power-law fit to PACS differential number counts in the 
form dN IdS °cS a . 



3. Resolved CIB fraction 

The added contribution of resolved galaxies provides a lower 
limit to the IR background and can be compared to di- 
rect measurements of the total CIB. After the discovery of 
the CIB (Pugetetal. 1996; Hauseretal. 1998), numerous au- 
thors attempted to directly measure its surface b r ightness from 
COBE/DIRBE maps (e.g. lLagache et ail l2000t iRenault et all 
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1200 It IWrightll2004l IDole et al.ll2006h Here we adopt the most 
recent revision of the Lagacheetal. (2000) DIRBE measure- 
ments, provided in Dole et al. (2006): 14.4 + 6.3, 12.0 + 6.9, 
and 12.3 ± 2.5 [nW m -2 sr _1 ], at 100, 140, and 240 fim, respec- 
tively. Interpolating between these values, one obtains a value of 
12.8 + 6.4 [nWrrT 2 sr" 1 ] at 160 //m. These direct measurements 
are still affected by large uncertainties, mainly due to difficulties 
in defining an absolute flux scale and in removing zodiacal light. 
An alternative estimate of the total CIB is obtained by integrat- 
ing the power-law extrapolation of our number counts, between 
0.01 and 1000 mJy. We adopt the slope derived on GOODS- 
N data for S < S break and an average between COSMOS and 
Lockman Hole at brighter fluxes (see Table [2] and Sect. [2J. The 
total CIB estimates from this extrapolation are ~12 and ~13 [nW 
irT 2 sr -1 ] in the green (100 fim) and red (160 fim) bands, within 
the errors of the direct measurements. 

The contribution to the CIB of individual GOODS-N sources 
detected above 3<x is vl v - 4.46 + 0.52 at 100 fim and 4.41 ± 0.62 
[nW or 2 sr^] at 160 fim, i.e., -31+4% and ~ 35 ± 5% of the 
direct DIRBE measurements (Fig. [2). Taking the effects of com- 
pleteness into account (see Sect. [2ji, we can compute the CIB 
fraction due to the general FIR-population, overriding source ex- 
traction losses. We limit the analysis to the 3cr detection limit 
again. The GOODS-N number counts cover the ranges 3-45 
mJy (at 100 fim) and 5.5 - 72 mJy (at 160 fim); the COSMOS 
field allows upper boundaries to be extended to 142 and 179 
mJy. Integrating the counts over the whole flux range, we obtain 
vl v = 6.36 + 1.67 and 6.58 ± 1.62 [nW m -2 sr 1 ], within 1.3cr 
and l.Ocr from the reference values at 100 and 160 fim, respec- 
tively (Fig.[2j». Th ese cor respond to - 45 ± 12% and ~ 52 + 13% 
of the IDole et a l. (2006) estimate. The derived uncertainty now 
also includes the effect of Poisson statistics and not only photo- 
metric errors. The bright end of number counts, covered only by 
COSMOS, gives a small contribution (~ 4 - 6%) to the total CIB 
surface brightness. 

The PACS detection rate of Spitzer 24 fim sources (S 24 > 20 
fiiy, Magnelli et al. 2009) is roughly 15%. It is possible to derive 
a deeper CIB estimate, through stacking on the PEP maps at 
the positions of all 24 ^m objects, including those not detected 
in the FIR (e.g. IDole et alH20 06). Uncertainties on the stacked 
fluxes are computed via a simple bootstrap procedure. The CIB 
surface brightness produced by 24 fim sources is 7.39 + 0.48 and 
9.57 + 0.71 [nWm" 2 s r -1 ] at 100 and 160 /xm, consistent with 
iBefhermin et al.1 d2010l) . and providing -51 % and -75% of th e 
total background, within 1 . lcr and 0.5cr from lDole et al.l (|2006). 



4. Discussion 

In the attempt to reproduce the observed ISO and Spitzer number 
counts, several authors built "backward" evolutionary models, 
including luminosity and/or density evolution, as well as differ- 
ent galaxy populations. In Fig.[T| we overlay recent models onto 
the observed PACS co unts. We include in this collection also 
the lLacev et al] d2009l) A-CDM semi-analytical model (SAM), 
complemented with radiative transfer dust repro cessing. 

The most successful models are the Valia nte et al.l ([2009), 
including luminosity-dependent distribution functions for 
the galaxy IR SEDs and their AGN contribution, and 
Rowan-Robinson (2009), employing analytic evolutionary func- 
tions without discontinuities and 4 galaxy populations. 

Gruppioni and Lag ache et al.1 (12004 ) overestimate the am- 
plitude of the number cou nts peak in both ban ds, while 
iFranceschini et al.l d2009l) and lLe Borgne et all (|2009) reproduce 
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Fig. 2. Contribution of GOODS-N galaxies to the CIB. Red 
empty triangles represent the sum of individually detected 
sources; black filled triangles come from number counts above 
the 3<r limit, and include completeness correction; blue upside- 
down empty triangles belong to stacking of 24 fim sources. Solid 
histograms denote the contribution of different redshift bins to 
the CIB, accounting for completeness correction. Hatched his- 
tograms are stacking results. Literature data include: DIRBE 
measurements revised by Dole et al. (2006, filled circles, lcr er- 
rors)' FIRAS spectrum (so lid lines above 200 ^m. lLagache et al.l 
ll999il2000l) ; lFixsen et alj d!998l) modified Black Body (shaded 
area). 



the counts fairly well only in one channel (100 and 160 fim, re- 
spectively). It is worth recalling that to date most of these models 
have been fine-tuned to reproduce mid-IR and sub-mm statistics, 
while a big gap in wavelength was affecting far-IR predictions. 
Most include a luminosity evolution oc (1 + z) ~ 3 , but the red- 
shift limit for this slope, the details of density evolution, or the 
adopted galaxy zoo vary significantly from author to author. 

Besides spanning a much wider range of observational data 
(UV, optical, near-IR luminosity functions, galaxy sizes, metal- 
licity, etc.), the SAM approach suffers here for a limited flexibil- 
ity in the choice of parameters, and significantly overestimates 
the bright end of PACS counts. Moreover, it cannot reproduce 
the peak, especially in the green band. 

Thanks to the rich ancillary dataset in GOODS-N, we split 
the far-IR number counts into redshift bins (Fig.[3}. This elabora- 
tion offers a remarkable chance to set detailed constraints on the 
evolution of the galaxy populations adopted in current recipes. 
This view highlights some new features and the main problems 
of the models under discussion. 

First, in the lowest redshift bin, 0.0 < z < 0.5, the differ- 
ential number counts, normalized to the Euclidean slope, mono- 
tonically increase as a function of flux, resembling the trend ex- 
pected for a non-evolving population of galaxies. The consis- 
tency of models with data weakens at longer wavelengths. At 
the bright end, above 10 mJy and 20 mJy at 160 fim, PEP cat- 
alogs include 43 and 18 objects in this redshift range. Based 
on template fitting (Rodighiero et al., 1201 Oh . half of thi s sam- 
ple shows a late-spiral SED (Sd-Sdm in the Polletta et al., 2007, 
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Fig. 3. Slicing the PEP GOODS-N population into redshift bins. Each PACS source has been associated to an optical/mid-IR coun- 
terpart as described in AppendixfA] See Fig.Q]and text for details on the models shown. 



template library) with typical luminosities above 10 11 L and a 
peak wav elength of ~130 urn. Nevertheless, at these luminosi- 
ties, |CharyjfcElbIl (HoOlb models, which adopt the local L - T 
relation, predict SEDs peaking around 90 yimi restframe, thus in- 
dicating that this class is indeed colder than typical local IR- 
luminous sources. This is reminiscent of the col d sources found 
by ISO FIRBACK and se rendipity surveys (e.g. Dennefel d et alJ 
l2005HStickel et alJl2007h . 

Significant differences between models arise at higher red- 
shift. The bin 1.0 < z < 2.0 shows the largest discrepancies: 
apparently, most of the overestimation of total number counts by 
some models originates here, suggesting that the adopted evolu- 
tion is either too protracted in time or too steep. 

We are now able for the first time to infer how much of the re- 
solved CIB peak was emitted at different cosmic epochs. Figure 
|2]depicts the amount of extragalactic IR background originating 
from the four redshift bins taken into account (see also Fig. [3}. At 
100 /mi, at the GOODS-N 3cr limit of ~3 mJy, roughly 80% of 
the resolved CIB is emitted by objects detected at z < 1, equally 
distributed below and above z = 0.5. The remainder belongs to 
more distant galaxies, with only 5% locked into individual ob- 
jects at z > 2. Most of the resolved 160 fim CIB is produced at 
higher redshift, with the 0.5 < z < 1.0 redshift bin dominating 
(-40%) the budget. 

At the PEP depth, the CIB resolved into individual sources is 
mainly due to luminous and ultraluminous IR galaxies (LIRGs 
10 11 < L IR < 10 12 L , ULIRGs L IR > 10 12 L ). The relative 
fraction of the two classes varies as a function of redshift, likely 
due to selection effects, with LIRGs providing 80-90% of the 
CIB resolved at 0.5 < z < 1 and ULIRGs dominating at higher 
redshift (60-70% at 1.0 < z < 2.0 and 100% above). 

Finally, Fig. [2] also reports the CIB contribution of 24 yum 
sources, as obtained through stacking (see Sect. [3]). The relative 
CIB fractions emitted in the four redshift bins change signifi- 
cantly at 100 yum, with respect to those obtained by integrating 
the observed counts. The redshift distribution of the CIB is now 
peaking between 0.5 < z < 1.0 and the relative contribution of 
z > 1 increases by more than 10%. On the other hand, at 160 
yum only the relative fraction of the highest redshift bin varies by 
more than 5%, while the others remain unchanged within 1 -2%. 



Also in the case of stacking, the 160 //m CIB relative redshift 
distribution is more populated at high redshift than the 100 //m 
one. As expected, this indicates that PACS galaxies with redder 
observed colors lie on average at higher redshift than bluer ones. 
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Appendix A: PACS data 

PEP science demonstration data include the GOODS-N and 
Abell 2218 fields. In addition to these, part of the other 
PEP blank fields have already been scheduled and observed: 
Lockman Hole and COSMOS (-85% of the planned depth). 
Observations of all fields were carried out by adopting the in- 
termediate speed (20 arcsec/s) scan-map mode. Table lATI lists 
the total exposure times already observed for these fieldqj. 

Data have been processed through the standard PACS reduc- 
tion pipeline, version 2.0.1328, within the HCSS environment 
(Ott et al. 2010). Additionally, we employed custom procedures 
aimed at removing of interference patterns, tracking anomalies, 
re-centering positional offsets, and mapping. 

Glitch removal i s based on multi-re s olutio n median trans- 
form, developed by IStarck & Murtaghl d!998l) to detect faint 
sources in ISOCAM data. The signal due to real sources and 
glitches show different signatures in the pixel timeline. These 
features are recognized using a multi-scale transform, separat- 
ing the various frequencies of the signal. Once the glitch com- 
ponents are identified, they are replaced by interpolated values 
in the pixel timeline. 

PACS photometers exhibit a noise with a roughly /" spec- 
trum at relevant frequencies. To remove the bulk of the noise 
we apply a "running-box" high-pass median filter to each pixel 
timeline, but mask the position of bright sources. The objects 
mask is produced iteratively during the reduction by detecting 
sources on the final map, and then we mask them in a double- 
pass mapping scheme. Testing shows that this masked filtering 
method modifies the fluxes of point-like source by less than 5%. 

Imperfections, drifts, and errors in the pointing accuracy of 
the Herschel satellite were corrected by re-centering the data 
on a grid of known 24 yt/m sources populating the fields. Such 
objects were stacked for all scan-legs in a given direction in a 
given map repetition (or for a subset of those, in the very large 
COSMOS maps). The stacking result was then used to compute 
the average offset to be applied to this set of scan-legs, for a given 
direction, in a given map repetition. This procedure also implic- 
itly corrects for small timing offsets between pointing informa- 
tion and data. Absolute, systematic astrometric offsets turned out 
to be as high as 5 arcsec, while relative corrections between in- 
dividual submaps are approximately 1 arcsec. 

Map reconstruction is done via simple image co-addition, 
based on a simplifie d version of the "drizzle" method 
(Fruchter & Hook 2002). Given the high data redundancy in the 
GOODS-N field, the drop size is set to 1/8 of the input ar- 
ray pixel size. This corresponds to 1/5 and ~ 1/4 of the output 
pixel size at 100 fim and 160 /im, respectively, thus reducing 
the correlated noise in the final map. Fields with lower redun- 
dancy were mapped by adopting a smaller drop size (1/4 of the 
input PACS array pixel size). Images produced from each ob- 
servation were weighted according to the effective exposure of 
each pixel and co-added to produce the final maps. The final er- 
ror map was computed as the standard deviation of the weighted 
mean. Owing to the nature of scan maps, correlations exist be- 
tween nearby pixels, in particular along the scan direction. These 
correlations are close to uniform across the final map, thus we 
derived a mean correlation correction factor which was then ac- 
counted for in the errors on the extracted fluxes. 



1 See the PEP web page for information about the other fields in the 
survey: http://www.mpe.mpg.de/ir/Research/PEP/ 

2 HCSS is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia 



PACS catalogs were extracted following two different ap- 
proaches, optimized for the different scientific aims of the PEP 
project. We perfor med a blind e xtraction using the Starfinder 
PSF-fitting code (iDiolaiti et al] l2000t) and a guided extrac- 
tion using 24 y um pr iors, following the method described in 
Magnelli et al. (120091) . The two methods provide similar results: 
fluxes extracted in the two cases are consistent with each other, 
the prior extraction leading to slightly deeper — although possi- 
bly biased — catalogs. The number counts presented in this pa- 
per were built on the blind catalog. Point spread function (PSF) 
profiles were extracted from the final science maps, and turn out 
to have an FWHM of -7.5 and -11 arcsec in the 100 /im and 160 
fim bands, respectively. Aperture corrections were characterized 
on calibration observations of the Vesta asteroid. Absolute flux 
calibration is based on yDra, cTau, and oCMa and makes use 
of the standard calibration file embedded in the PACS pipeline 
(version 2.0.1328). Typical absolute flux calibration errors are 
~ 10% and include uncertainties on instrumental characteriza- 
tion, PSFs, and reference stars analysis. 

Noise in PACS maps was measured with random aperture ex- 
tractions on residual images and compared to the observed S/N 
ratio for the detected sources. The r.m.s. values thus obtained 
include both instrumental noise and confusion noise due to un- 
detected sources (i.e., below the 3cr threshold). This measured 
lcr noise is 1.00 mJy at 100 //m and 1.90 mJy at 160 fim for 
GOODS-N. Table IA. II includes the noise properties of the SDP 
fields, as well as the Lockman Hole and COSMOS. 

To quantify the reliability of extracted fluxes, the level of in- 
completeness and the fraction of spurious sources, Monte Carlo 
simulations were performed, creating 500 images and adding 20 
artificial objects onto science maps each, for a total of 10000 
sources. Input and output fluxes are consistent with each other 
within a few percent. Completeness is defined as the fraction 
of sources that have been detected wi th a photometric accuracy 
of at least 50% dPapovich et al.ll200 4). Spurious sources are de- 
fined as those extracted above 3cr with an input flux lower than 
3cr{Image). The latter is consistent with the spurious fraction in- 
ferred by blindly extracting from inverted maps. The GOODS-N 
blind catalog reaches 80% completeness at -5.5 mJy and -11.0 
mJy in the two bands, and a 30% fraction of spurious detections 
at -2.5 and -7.0 mJy, in green and red respectively (see Table 
lATI) . 

The GOODS-N field benefits f rom an exten s ive multi- 
wavelength coverage. Adopting the iGrazian et al] ([2006) ap- 
proach, the PEP Team built a reliable multi-wavelength, PSF - 
matched database, including ACS bviz dGiavalisco et alj |2004). 
Fla mingos JHf$ and Sp itzer IR AC data. Moreo ver, MIPS 24 
Hm (Magnelli et al. 2009) and Ba rgeretal.1 d2008h deep U, Ks, 
and spectroscopic redshifts have been added. When no spectro- 
scopic redshifts were available , photometric redshift s have been 
derived using the EAZY code (Bram mer et alj [2008). 

The 100 and 1 60 jum blind catalogs were finally linked to this 
multi-wavelengt h catalog through a three-band maximum likeli- 
hood procedure (Sutherland & Saunders 1992), starting from the 
longest wavelength available (160jum, PACS) and progressively 
matching 100 ^m (PACS) and 24jum (Spitzer/MIPS) data. Table 
[AT] summarizes the main properties of blind catalogs, as well as 
the results of the maximum-likelihood match in GOODS-N. 

Confusion is a major concern in deep wide-beam obser- 
vations, like far-IR or sub-mm imaging of blank fields. PEP 
Herschel observations are not dispensed from confusion: the 
high density of detected sources hinders the extraction of fainter 
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objects, in the s o-call ed "source de nsity confusion crite rion" 
(SDC, Dole et al. l2003h . Adopting the lLagache etail (120031) def- 
inition of beam (i.e., D. = 1.14 x Q 2 FWHM ), the source density 
in PACS blind catalogs is 40 beams/source at 100 fim and 18 
beams/source at 160 /mi. This indicates that PEP GOODS-N is 
already hitting the SDC limi t in the PACS red band, estimated 
to be 16.7 beams/source bv iDole et alJ (120031) . while the green 
channel is not affected. By extrapolating integral number counts, 
we estimate that this limit will be reached at ~2.0 mJy at 100 
yum and ~4.7 mJy at 160 /vm. Nevertheless, this limit is known 
to be rather conservative, and some authors have already shown 
that source extraction can be reliably carried out to levels as low 
as 10 beam/source under favorable conditions (e.g. at 24 //m, 
iMaenelli et afll2009h . 
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Notes. {a) The PEP Lockman Hole field is the ROSAT-HRI and XMM field. (b) COSMOS will reach 213 hours of integration at full depth. 
Table A.l. PEP fields: total exposure times, noise properties, flux levels for 80% completeness and 30% spurious fraction, statistics 
of blind catalogs, and results of maximum-likelihood match to the multi-wavelength ancillary catalogs (labeled "multi", GOODS-N 
only). 
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Notes. Fluxes are provided in [mJy]. For each field/wavelength, we list counts in units of [deg 
include both Poisson statistics and propagation of photometric uncertainties. 

Table 1. PEP number counts, normalized to the Euclidean slope. 



mJy 1 5 ]. Errors are given as relative fractions, and 



